#include "udf.h"

DEFINE_PROFILE(pressure_profile,t,i)
{
  real x[ND_ND];
  real y;
  face_t f;
  begin_f_loop(f,t)
  {
     F_CENTROID(x,f,t);
     y = x[1];
     F_PROFILE(f,t,i)=0.5*y;
  }
  end_f_loop(f,t)
}